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ABSTRACT 

We present the first spectral energy distributions produced self-consistently by 2.5D 
general relativistic magneto-hydrodynamical (GRMHD) numerical simulations, where 
radiative cooling is included in the dynamical calculation. As a case study, we focus on 
the accretion flow around the supermassive black hole in the Galactic Centre, Sagit- 
tarius A* (Sgr A*), which has the best constrained physical parameters. We compare 
the simulated spectra to the observational data of Sgr A* and explore the parameter 
space of our model to determine the effect of changing the initial magnetic field config- 
uration, ion to electron temperature ratio Ti/T e and the target accretion rate. We find 
the best description of the data for a mass accretion rate of ~ lO -9 M0/yr, and rapid 
spin (0.7 < a* < 0.9). The submillimcter peak flux seems largely independent of initial 
conditions, while the higher energies can be very sensitive to the initial magnetic field 
configuration. Finally, we also discuss flaring features observed in some simulations, 
that may be due to artifacts of the 2D configuration. 

Key words: accretion, accretion discs - black hole physics - MHD - radiation 
mechanisms: thermal - methods: numerical - Galaxy: center. 



1 INTRODUCTION 

The best studied low-luminosity active galactic nucleus 
(LLAGN) is the supermassive black hole in the centre 
of the Milky Way, Sgr A*, disco vered originally via it s 
strong radio continuum emission bv lBalick fc Brownl l|l974l) . 
Sgr A* is a unique system because of its proximity com- 
pared to the centres of other galaxies, and a multitude 
of intensive single- and multi-wavelength campaigns hav- 
ing been conducted from the radio through high-energy 7- 
rays over t he last decades (see references in Me lia fc Falckel 
l200ll and iGenzel et all l2010l ). These studies provide re- 
markably stringent constraints on Sgr A*'s properties. The 
current best mass, distance and mass accretion rate val- 
ues are M = 4.3 ± 0.5 x 1O 6 M , D = 8.3 ± 0.4 kpc 
and 2 x 10~ 9 M R yr" 1 < M < 2 x 10~ 7 M R yr~\ respe c- 



tively jReidl|l993l:lGhez et al.ll2008l:lGillessen et al.ll2009al lbl; 
iBower et al.ll2005l : iMarrone et al.ll2007h . 

The above constraints make Sgr A* the perfect candi- 
date to test theoretical models of the accretion processes at 
low accretion rates. In particular, high precision data from 
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close to the event horizon allow us to study the detailed 
physics of the accretion flow and potential jet launching, 
in the extreme regime where gravity is important. Further- 
more, because Sgr A* is representative of the majority of 
SMBHs today, often lurking below the detection threshold 
of even our most sensitive telescopes, we can use it to get 
a handle on the contribution from these very weak AGN to 
their host galaxies. 

In the attempt to understand the nature of accretion 
flows, many semi-analytical models have been developed. 
Accretion disc models li ke the thin disc developed by 
IShakura fc Sunvaevl ll 19731). advection-dominated accretion 
flow f ADAF; iNaravan fc Yil Il994l ) advection-dominated 
inflow -outflow solutions (ADIOS; iBlandford fc Begelmanl 
1999) , convection-dominate d accretion flow ( CDAF; 
Quataert fc Gruzinovl [2000b) or Bondi accretion dMelial 
19921 . 1994) have been successful at fitting data from many 
sources. However, by nature of the semi-analytic approach 
which cannot model turbulence, none of these models can 
accurately address the role played by the magnetic field in 
the dynamics of accretion disc via the magnetorotational 
instability fMRI; iBalbus fc Hawlevlll99ll ). 

Recent breakthroughs in parallel numerical simulations 
now allow extensive fluid simulations with full general rela- 
tivity (GR) and magneto-hydrodynamic (MHD) treatment 
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in reasonable timescales. Such simulations have finally pro- 
vided the opportunity to study the dynamical properties of 
an entire system in a complementary manner to the semi- 
analytical models. In particular, various GRMHD codes 
have been employed by several groups over the past few 
years, in order to perform detailed th eoretical studies of 
the accretion flow around Sgr A* (e.g . [ Dexter et al.l 20091; 



Moscibrodzka et al. 2009 



2010: Shchcrbakov et al 



.Hilburn et al.l 
20101; iMoscibrodzka et alj l201ll ; 



Dexter et al|l2010l 



Dexter fc Fragilell2012l ; lDolence et alj|2012h . with the goal of 



reproducing Sgr A*'s spectra. These studies all share a com- 
mon approach in which the radiative losses are not included 
in the simulations themselves, but rather first a dynamical 
model is calculated in GRMHD, and then the final outputs 
are fed into a separate post-processing routine to calculate 
the resultant spectrum. These studies justified ignoring the 
inclusion of cooling because Sgr A* is so underluminous that 
radiative losses are likely not strong enough to affect the dy- 
namics of the system. 

In a companion paper l|Dibi et al.ll20l3 . hereafter iDllh . 
we assess, for the first time, the importance of the radia- 
tive cooling in numerical simulations of Sgr A* by using 
Cosmos++, an astronomical fluid dynamics code that takes 
into account radiative losses self-consistentl y in the dynam- 
ics l|Anninos et al ] |2005l : iFragile et alll2012h . We show that, 
for Sgr A*, cooling effects on dynamics can indeed be ne- 
glected. However, the effects of cooling at higher accretion 
rates (relevant for most nearby LLAGN) are not negligible. 

In this paper, we describe the implementation and re- 
sults from the cooling routines used in the simulations of Sgr 
A* presented in|D12j, and present the first self-consistently 
calculated spectra in order to explore the new parameter 
space. We examine the influence the spin and the initial mag- 
netic field configuration have on the simulated spectra, and 
compare to the previous non-cooled calculations. Although 
we find that self-consistent treatment of radiative losses is 
not important for the case of Sgr A*, we demonstrate that 
it will be for most nearby LLAGN. 

Section [2] describes the observational constraints on 
Sgr A*. In Section [3] we present the derivation of emissiv- 
ity expressions that generate the cooling rates employed in 
Cosmos++. We also discuss numerical limitations of our sim- 
ulations, as well as the assumptions made when generating 
the spectra. In Section [4j we present the results obtained 
by comparing our spectra to observational data, and in Sec- 
tion [S] we discuss them. Finally, in Section [U we summarise 
our conclusions and suggest future improvements. 



2 OBSERVATIONAL CONSTRAINTS 

The radio spect r um of S gr A* (|Serabvn et al.l 1 19971 ; 



iFalcke et all 1 19981 ; I An et alj|2005h sh ows a slight chang e 
in the spectral index above 10 GHz ( Falcke et al. 1998), 
peakin g in a so-called submtllimeter bump. Aitken et al.l 
(2000) reported the detection of linear polarisation from 
this submm bump, that is not observed at longer wave- 
lengths. This change implies that the radio emission and the 
submm bump originate from distinct but contiguous regions 
in the system. Very long baseline interferometry measure- 
ments have limited the size o f the submm-emitting region 
to be ~ 4 Schwarzschild radii |Doeleman et alj|2008l ). 



Limited angular resolution and sensitivity in the far- 
and mid-infrared bands makes it impossible to distin- 
guish between the emission from Sgr A* and the sur- 
rounding sources, resulting in no observations in the far- 
infrared band and only upper-l imits on the flux in the 



infrared (e.g 


. iDavidson et al.l 1 19921; iHerbst et al. 


19931; 


Stolovv et al. 


1996;lTelesco et al.lll996l; iMenten et al. 


1997; 


Schodel et al. 


2007. 20111) and in the X-rav (Baganoff et al.1 


2003: Bclangcr et al. 20061). show a quiescent state lower 



in flux than the submm bump, where the radiation power 
peaks. However, a few times a day, Sgr A* experiences 



rapid increases in t he ne a r-infrared flux ([Hornstein et al 



20021 ; iGenzel et all 120031; iGhez et al l l20oi 



2006 



2011 



Eckart et al 



20081; lYusef-Zadeh et al.l l200sl ; iDodds-Eden et al 
Haubois et al.l 120121') . whe re brighter flares (> 
Dodds-Eden et al.1 l201ll ) are often associated 
X -ray flares 



10 mJy 

with simultaneous 



Goldwurm et al l 120031; [Porquet et all 2003; iBelanger et al.l 



jBaganoff et al.1 l-'oOl: 
I. 



20051 ; IPorquet et al.ll200sir 

These multi-wavelength observatio ns provide tight con- 
strain ts on the physics of Sgr A*. Sha kura fc Sunvaevl 
|l973t ) type radiatively efficient, thin disc models are ex- 
cluded as they predict an observed infrared flux several or- 
ders of magnitude hi gher than the upper-limits obtained 
jFalcke fc Melialll997h . The presence of linear polarization 
and the constraints on the Faraday rotation currently limit 
the mass accretion rate of Sgr A* to be much smaller 
than the Bondi accretion rate, of the order of ~ 10~ 8 M (?) 
|Aitken et al.ll2000l ; iBower et al"1l2003l ; iMarrone et~aI1l2007l ) . 
Such a low accretion rate in fact exclude s the "class ical" 
ADAF (jNaravan et al.1 Il998l ) and Bondi (|Melial 1 19921 ) ac- 
cretion models for Sgr A *, as these models invoke higher 
accret ion rates (see, e.g.. lAgoll2000l ; IQuataert fc Gruzinovl 
l2000al ). In the meantime, many other models have been 
developed that are still consistent with the current lim- 
its. Radiatively ineffic i ent accretion flow models (RIAF ; 
i Blandford fc Begehnar] Il999l ; IQuataert fc Gruzinovl boOObl ; 
lYuan et al.l 120031 ) argue that the submm emission is pro- 
duced via synchrotron radiation from a thermal distribu- 
tion of electrons, in the innermost region of the accretion 
flow, whic h could also be synonymous with the b ase of 
the jets l|Falcke fc MarkofJ [2OO0I ; I Yuan et al.1 12002| ). This 
synchrotron emission is then inverse Compton upscattered 
by these same electrons, resulting in a second peak that 
contributes to the X-ray emission during flares. The ra- 
dio emission can originate from eith er a non-thermal tail 
of electrons produced in a RIAF (|Yuan et al.l 120031 ) or 
from predo minantly thermal electr ons within a mildly rela- 
tivistic jet (IFalcke fc Mar kofl 2000[). B ased on o b servat ions 
with Chandra (|Baganoff et al.1 120031 ). IQuataert! (|2002l ) ar- 
gues that the faint quiescent X-ray emission is from thermal 
bremsstrahlung, originating in the outer region of the accre- 
tion disc. 

While successful at producing a general description of 
the data, all the above semi-analytical models lack a self- 
consistent MHD description of the accretion flows. Although 
they invoke a viscosity to account for the outward angular 
momentum transport, in the accretion discs, they do not 
explicitly calculate it, nor do they account for the presence 
of magnetorotational in stability driven accretion processes 
ijBalbus fc HawleyHl99ll ). GRMHD simulations are thus an 
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ideal framework to examine the nature of accretion flows 
around black holes, and to test the above scenarios for Sgr 
A*'s emission in particular. 

The GRMHD simulations presented in this work model 
only the innermost part of the accretion flow around Sgr A*, 
where the submm bump is produced. Therefore the radio 
(jets or outer accretion inflow) and X-ray (outer regions of 
the accretion disc) emission cannot be fitted by this present 
work and we focus our results on fitting the submm emission 
of Sgr A* in the quiescent state. However, we use the IR/X- 
ray emission as upper limits to define the feasibility of our 
fits. 



3 METHODS 

The general setup of our simulations i s similar to that of 
other groups llMoscibrodzka et aDl 2009: Hil burn et al.ll20icl : 
iMoscibrodzka et al.ll201lh in order to facilitate comparison. 
The simulations start with an initial torus of gas, seeded 
with a magnetic field, around a compact object situated at 
the origin. The mass of the central object is set to the mass 
of Sgr A* (Mbh = 4.3 x 10 6 M Q ) and the initial density 
profile inside the torus is chosen to produce the target mass 
accretion rate at the inner grid boundary. We let the sim- 
ulation evolve until inflow equilibrium is established in the 
inner disc, where a jet, along the rotation axis, and an out- 
ward flowing wind, which we call t he cor ona, over and under 
the accretion disc, are formed. See|Pl3 for more details. 

To generate spectra, we had to ensure that the radia- 
tive emissivities used were physically accurate and, when 
integrated, produce the same cooling functions as origi- 
nally included in Cosmos++. We thus consider the following 
cooling processes: bremsstrahlung, synchrotron and inverse- 
Compton. Since we are investigating physical processes oc- 
curring very close to a black hole, we also must include spe- 
cial and general relativistic effects on the radiative emission. 
In the following, we describe the adopted emissivity expres- 



3.1 Radiative cooling 

Cosmos++ uses as a cooling functi on the following total cool- 
ing rate for an optic ally thin gas ^Fragile fc M eier 2009], also 
see lEsin et ai1ll996l ): 



A = ribr,c q br + Vs,c q s 



(1) 



where q^ r and gj" are respectively the bremsstrahlung and 
synchrotron cooling terms and 7]b r ,c and r)s,c are Compton 
enhancement factors. These rj factors are modified exponen - 
tial function of the Compton parameter y IjEsin et all 1996). 

The brems strahlung cooling rate is taken from 
lEsin et all (1 19961 ) (equations (7) to (9) in this paper) 



Qbr = lei + lee + <?± 



(2) 



where q~ t , q^ e and q± represent the cooling due respec- 
tively to electron-ion and positron-ion, electron-electron and 
positron-positron, and electron-positron processes. 

The synchrotron cooling rate is a sum of optically thick 
and thin emission (equation (14) in lEsin et alj (| 19961 )) 



<ls = 



2nkT 
He 2 



v dv + 



,(u) di 



(3) 



where k is the Boltzmann constant, T is the temperature of 
the electrons, H is the local temperature scale height, c is the 
speed of light, v c is the critical frequency at which the opti- 
cally thick and thin emissivities are equal a nd e a (f) is the to- 
tal a ngle- averaged synchrotron emissivity (|Fragile fc Meier] 
l2009h . 

The cooling rates are important to evaluate the radia- 
tive losses at each time step of the simulation. Whereas the 
emissivities, from which these cooling rates are derived, are 
the critical quantities we need to produce the spectra. 



3.1.1 Bremsstrahlung 

Similar to lEsin et al.1 l|l996l ). we used IStepnev fc Guilberd 
i| 19831 ) expression for the thermal relativistic electron-ion 
bremsstrahlung emissivity 



dE e 



dVdtdw 



aw 



(4) 



where w = hv/m e c 2 is the dimensionless photon energy, 
N e (rf) = N e j 2 /3 exp(-j/6)/6K 2 (l/6) is the Maxwellian- 
Juttner electron energy distribution, 8 — kT/m e c 2 is the di- 
mensionless temperature, and K2 is a modified Bessel func- 
tion. N e and N p are the electron and proton number densi- 
ties, respectively. 

The electron-ion cros s-section needed in equation if!"} 
can be expressed following iBlumenthal fc Gouldl (|l970l ) 



da 16 2 2 * w 3 f w 
w - — = — Z ar 1 h - — 

7 4 V7 



dw 



hi 



2 (7 2 — 7«j) 



(5) 



where Z is the ion's atomic number, a is the fine structure 
constant and ro = 2.8179 x 10~ 13 cm is the Compton radius. 

Combining equations (|4]) and ((5| gives the following ex- 
pression for electron-ion bremsstrahlung emissivity 

dE ep 16 2 2 f °° a 2 I 2 . 3 2 

—— = N e N p c Y Zar J^J (7 - JW + 

/ 2( 7 2 -yw) A e~~<l 9 , , s 
x In — !— £ - - -pr-drt (6 

Integrating equation ((6} over frequencies leads to the 
bremsstrahlung cooling rate described by equation ©. 

When compared to electron-ion bremsstrahlung interac- 
tion, electron-electron and electron-positron bremsstrahlung 
are negligible. Therefore their contribution to spectra has 
been ignored in the present work. Moreover, in the region 
of interest in our numerical simulations (i.e. r < 15 r g , see 
Section 14. ip , bremsstrahlung as well as Comptonization of 
bremsstrahlung have a smaller contribution to the overall 
emission in comparison with both the synchrotron and the 
synchrotron self-Compton processes. However, although in 
the case of Sgr A* bremsstrahlung emission is negligible, it 
will have an important contribution to the spectra of other 
LLAGN. 



3.1.2 Synchrotron and synchrotron self-Compton 

The total angle-averaged opt ically thin and thick syn- 
chrotron emissivities given bv iFragile fc MeTel] (fioogh arc 
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only valid within a certain range of temperatures. There- 
fore, rather than using them to account for the synchrotron 
contribution to the spectrum, we decided to start from first 
principles to express a more g eneral expression of t h e syn - 
chr otron emissivity. Foll owing iRvbicki fc Lightmanl (1 19861 ) 
and lde Kool et alTl) 1989^ 1 . we have: 



(1 - exp(-p„i?)) 



(7) 



where r\ v is the emission coefficient, /i„ is the absorption 
coefficient and R is the size of the homogeneous emitting 
volume. 

K nowing the synchrotron radiat ion field 7 S (y) and fol- 
lowing IChiaberge fc Ghisellinil (|l999l ) , the synchrotron self- 
Compton emissivity, in units of erg/cm 3 /s/Hz/st can be ex- 
pressed as 



e c {vi) = — 



dvg 



"' 7 iv( 7 )/K*i)-/*M 



7 2 £ 



(8) 



where is the frequency of the incident photons, v\ is the 
frequency of the scattered photons, f(uo, v\) is the spectrum 
produced by the single electron, scattering monochromatic 
photons of frequency u , (3 = 2, v™" 1 and v™ ax are the 
extreme frequencies of the synchrotron spectrum, and 71 
and 72 are 



71 



(- 

\4voJ 



1/2 



1 ^min 



72 



^fmax , 



3 m e c 



4 hvo 



(9) 



(10) 



The mean free path of inverse Compton scattering being 
larger than the simulation region in all our cases, we ignore 
multiple scatterings along the line of sight. 

For the synchrotron self-Compton radiation field from 
a homogeneous volume of size R, equation ([8)l leads to the 
following emissivity expression: 



Ic(vi) = e c {vi)R 



(11) 



In the framework of the simulation, R represents the size of 
a zone. 

Integrating, over frequencies, the synchrotron and syn- 
chrotron self-Compton radiation fields, I 3 (y) and I c {v), 
leads to the synchrotron cooling rate, with the Compton 
enhancement factor r} s ,c Q7 ■ 

These expressions integrate to exactly the formulae used 
for the cooling rates within Cosmos++. 



3.2 General Relativistic Radiative Transfer 

The synchrotro n emission and absorp tion coefficients (equa- 
tions (4)- (5) in Ide Kool et al] 1 19891 ). and the synchrotron 
self-Compton emission coefficient (equation [8]) describe the 
emitted spectrum from any zone in the simulation. A ra- 
diative transfer calculation is necessary to transform this 
into the spectrum as seen by a distant observer. Due to 
both strong gravitational lensing and redshifts, and Doppler 
beaming in the vicinity of the black hole where most of the 
luminosity is produced, this calculation must be done in full 
GR. 



The GR calculation is done using ray tracing. Start- 
ing from a distant observer's hypothetical detector, rays are 
traced backwards in time toward the black hole assuming 
they are null geodesies (geometric optics approximation) , 
using the public code geokerr described in iDexter fc Agoll 
l|2009h . In the region where rays intersect the accretion flow, 
the radiative tra nsfer equation is solv ed along the g e odesic 
|Broderickll2006r ) in the form given in iFuerst fc Wui (j2004T l 
using the code grtrans (|Dexterl201ll '). which then represents 
a pixel of the image. This procedure is repeated for many 
rays to produce an image, and at many observed frequencies 
to calculate the spectrum. 

Both gravitational redshifts and Doppler shifts lead 
to differences between observed and emitted frequencies. 
Emission and absorption coefficients are then interpolated 
both spatially between neighbouring zones to points on the 
geodesic, but also logarithmically in frequency to the emit- 
ted frequency corresponding to the desired observed fre- 
quency. 



3.3 Assumptions and numerical limitations 

All of our models assume a thermal plasma. This plasma is 
described with a Maxwell- Jiittner energy distribution with 
temperature T e , characterised by a fixed fraction of the ion 
temperature Ti. This approach is standard for MHD simula- 
tions since it would be difficult computationally to simulate 
two interacting plasmas. The ion temperature is calculated 
via the ideal gas law. Since the internal energy of the plasma 
is dominated by the ions, the cooling function used in the 
simulation is that of the ions. The assumption is made that 
the temperature of the electrons, T e , needed in the calcu- 
lation of the cooling rate (since the cooling processes that 
we are considering all involve electrons), is simply related 
to Ti by a fixed factor. To get T e , we assume that some 
process is coupling the two temperatures. In the case where 
the ratio is 1, we assume that the two temperatures are 
coupled via a perfect process. When Ti > T e this process is 
assumed to be imperfect. There is no reason why the plasma 
remains at a fixed temperature ratio throughout its evolu- 
tion. However, studies have shown that allowing this ratio to 
be space- and time-dependent do not dramatically change 
the resulting simulations (based on unpublished work by 
Dexter). Therefore the assumption of a fixed Ti/T e may be 
a reasonable approximation. 

We also assume that the radiation escapes freely from 
the system. The whole system is optically thin to syn- 
chrotron self-Compton emission while, for the calculation of 
the synchrotron, we consider the appropriate optical depth 
of the gas at a given location and time, which depends on the 
state of the plasma. This approximation takes into account 
the optical depth without performing radiative transfer, and 
is valid as long as the (assumed thermal) peak of the radi- 
ating particle distribution corresponds to energies greater 
than the self-absorption frequency, which is almost always 
the case for the regions under study. 

Numerical caps and floors are a necessary limitation of 
most MHD simulations to prevent the codes from crashing 
in regimes where the values are too large/small. In our sim- 
ulations, floors have been applied on the matter and energy 
density values a zone is allowed to reach. In addition, a cap 
has been imposed on the magnetization of the fluid, as mea- 
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sured by Pb /(p+pe), where Pb is the magnetic pressure, p is 
the density and e is the internal energy. Similar floors must 
also be applied when post-processing the simulations. Re- 
gions of simulation that are under-dense and under-energetic 
(for instance inside the jets), as well as strongly magnetically 
dominated zones (for instance near the black hole) are not 
included in the spectral calculations, as they have reached 
these numerical floors. 



3.4 Magnetic field configuration 

Magnetorotational instabilities in MHD simulations are 
driven by weak poloidal magnetic field loops seeded in the 
initial torus. However not much is known about the magnetic 
field configuration in accretion discs around black holes. In 
the case of Sgr A*, most groups model it with one loop across 
the initial torus. In our work, we tested the results using two 
different configurations: a single set of poloidal loops (here- 
after the 1-loop model) centred on the pressure maximum 
of the torus and following contours of pressure/density; we 
also run simulations with four sets of poloidal loops (here- 
after the 4-loop model) spaced radially, with alternati ng fie ld 
directions in each successive loop (see Section 5.1 of ID 12T) . 

Because of the stochastic nature of MRI-generated tur- 
bulence and magnetic reconnection, MHD simulations of ac- 
cretion discs can show significant variability. In addition, our 
axisymmetric simulations show violent flaring events trig- 
gered by reconnection. The impact of such events on the 
emission can be extreme, especially for a few very brief X- 
ray flares f see I5.3f l . 



3.5 Spectral Energy Distribution 

We let our simulations run for 7 orbits, where we refer to the 
circular orbital period at r = r centre or t or i, = 1.67 x 10 4 s. 
The simulations reach their targeting mass accretion rates, 
after their peak value and before returni ng to their back- 
ground rates, between 2.5 and 3.5 orbits (|D12r ). To repro- 
duce the quiescent state of Sgr A*, we take the median value 
of the 50 individual spectra in this interval. We do not use 
a time averaged SED, in order to not overweight the likely 
unphysical flaring episodes. Figure \T\ shows that a simple 
time-average gives too much weight to the flaring events, 
increasing the flux up to an order of magnitude in the X- 
ray, compared to the median. The shadow region is the "1- 
sigma" variation about the median. It represents the limits 
within which 68% of the spectra fall. For each 50 individual 
spectra, the eight highest and the eight lowest data points 
in each spectra energy bin have been dropped. 

FigurefSJpresents a sample simulated broadband spectra 
of Sgr A*. The first bump from the submm band to the near- 
infrared band is due to thermal synchrotron radiation while 
the second bump in the X-ray is from upscattered submm 
seed photons via Inverse-Compton process. 



3.6 Parameter-space 

Each model used to simulate Sgr A* is described in terms 
of the following five parameters: the configuration of the 
magnetic field B, the spin of the black hole a*, the ion-to- 
electron temperature ratio Ti /T e , the mass accretion rate M 




10° 10 9 10 10 10" 10 12 10 13 10 14 10 15 10 16 10 17 10 18 10 19 10 2 ° 
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Figure 1. Broadband spectra of the reference simulation 
B4S9T3M9C computed from all time steps in the interval 2.5- 
3.5 orbits. Time-averaged values (blue, dash-dash) are compared 
with median values (black, solid). The flaring events we see in 
our simulation have too much weight in the time-averaged val- 
ues, which is why we have chosen to use median values to repre- 
sent typical flux densities. Shadows represents the 1-sigma vari- 
ability of the simulated data (see Section l3.5|l . Observation data 
of Sgr A* (pink ) show average quiescent s pectrum publi s hed in 
iMelia k. Falckej l|200ll l. submm data f r om iMufioz et al.l l|2012l l 
and mean infrared from ISchodel et al. The X-ray is an 

average quiescent flux from iBaganoffet al] l|2003h . 
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Figure 2. Broadband spectra of the B4S9T3M9C simulation 
presenting the synchrotron (orange, dash-dot-dot) and the syn- 
chrotron self-Compton (blue, dash-dash) components of the radi- 
ation and the resulting total emission (black, solid). 



and enabling (C) or disabling the cooling function. A sixth 
parameter, the inclination angle i at which the system is 
viewed from Earth is used in the ray-tracing program and 
has been also studied. Low inclination angle corresponds to a 
face-on situation while high one is edge-on. Table [1] presents 
our parameter-space. 

As explained in Section [2] there are tight observational 
constraints on the mass accretion rate from linear polarisa- 
tion measurements in the submm band. We therefore im- 
pose boundaries of 2 x lO _9 M yr _1 and 2 x lO~ 7 M yr _1 
to our corresponding parameter. On the other hand, despite 
all the data that have been gathered on Sgr A* over the 
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Table 1. Free parameters explored in the simulations 



Parameter 


Values 


B loops (N) 


1, 4 


Spin (a,) 


-0.9, 0, 0.5, 0.7, 0.9, 0.98 


Ti/T e 


1, 3, 10 


Target M[ M Q /yr] 


10" 9 , 10" 8 , 10" 7 


Cooling function 


ON (C), OFF 


Inclination angle i [deg] 


5°, 45°, 85° 



years, its spin is still an unknown parameter. Therefore we 
have decided to explore a wide range of possible spin values, 
from a non-spinning black hole case (a* = 0) to a maxi- 
mum spinning (a* = 0.98), as well as a retrograde-spinning 
(a, — —0.9) one. Finally we allow the ion-to-electron tem- 
perature ratio vary between 1 (efficient coupling between 
electrons and ions) and 10 (relatively inefficient coupling), 
which is the same range explored in earlier works. 

Throughout this pape r, we follow the same naming 
convention defined in |P12l to designate simulations. Each 
simulation name refers to the parameter-space explored in 
that model. For example, B4S9T3M9C means that the ini- 
tial torus is seeded with a 4-loop poloidal magnetic field, 
the spin of the black hole is set to a, = 0.9, the ion-to- 
electron temperature ratio to 3, the mass accretion rate to 
M — 1O _9 M0 yr _1 and the cooling function is enabled. Ta- 
ble [2] presents a overview of all of our simulations. 

B4S9T3M9C at an inclination angle of 85° is the closest 
solution to those found by previous works when attempting 
to fit Sgr A*'s data with simulated SEDs, except for the 
enabling of the cooling function, and the initial magnetic 
field configuration, fn the following we have chosen this set 
of initial parameters to be our reference simulation, from 
which we have explored our parameter space. 



4 RESULTS 

4.1 Geometry 

To compute our spectra, we had to ensure that only regions 
of the simulations that have reached inflow equilibrium are 
contributing to the emission. The reason is that properties of 
the parts of the simulations that have not reached the inflow 
equilibrium are strongly dependent on the arbitrary, initial 
conditions. As a consequence, only radiation from a region 
lying between t he eve nt horizon and 15 r g are accounted for 
m our spectra (|D12|) . Figd] illustrates this selection. This 
geometrical restriction is consistent with our aim to only fit 
the submm bump, which is believed to originate very close to 
the black hole. Indeed, Figure 0] shows that this inner region 
of the accretion flow accounts for the bulk of radiation in 
our simulation. 

4.2 Exploring the parameter space 

For the first time, we are able to generate consistent spec- 
tra from GRMHD simulations that can be compared to ob- 
servations in a robust way, as no post-processing scaling is 
possible when the cooling function is enabled. To assess the 
importance of self-consistent treatment of radiative losses 



Mass density 




m to so so IOC 

Radius ffljj 

Figure 3. Delimitation of the different regions of our simulations 
considered in the simulated spectra. The limits of each region are: 
inner region [event horizon < r < 15 r g ] - disc [15 r g < r < 70 r g 
: 0.36tt < < (1 - 0.36)tt] - outer region [70 r g < r < 120 r s : 
0.36tt < 6 < (1 - 0.36)tt] - corona [15 r g < r < 120 r g : 0.08tt < 
9 < 0.36vr] - jets [15 r g < r < 120 r g : 0.004tt < 6 < 0.08tt]. The 
corona and jet regions are symmetric with respect to 8 = n/2. 
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Figure 4. Broadband spectra comparing the emission from every 
part of the simulation as shown in Figure [3] For clarity, the 1- 
sigma variability has been omitted. 



on the resulting emission, we have compared SEDs from 
simulations with the same set of initial parameters (4-loop 
model, a* = 0.9 and Ti/T e = 3), with and without radiative 
cooling. This comparison was done for three different target 
mass accretion rates. Figure O presents the six spectra ob- 
tained. When enabling the cooling function, there is a clear 
trend of increasing importance of the effect with increasing 
mass accretion rate. While at a mass accretion rate of 10 -9 
Mq /yr, both spectra of cooling and non-cooling simulations 
are similar, significant differences (up to two orders of mag- 
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Table 2. Description of simulation parameters 



Simulation B loops (TV) Spin (a,) Ti/T e Target M statistical av. M Cooling Resolution 



B4S9T3M9C 


1 


0.9 


3 


10~ 9 


2.36±1.54 x 10 9 


ON 


256 x 256 


B4S9T3M9Ct^ 
B4S9T3M9Ct£ 


1 


0.9 


3 


10 -9 


3.05 ± 2.22 x 10~ 9 


ON 


256 x 256 


4 


0.9 


3 


10" 9 


1.95 ± 2.62 x 10~ 9 


ON 


256 x 256 


B4S9T3M9CeJ 


4 


0.9 


3 


10 -9 


3.00 ±2.00 x 10- 9 


ON 


512 x 256 


B4S9 


4 


0.9 


- 


- 


- 


OFF 


256 x 256 


B1S9T3M9C 


1 


0.9 


3 


10" 9 


7.93±9.27 x 10~ 9 


ON 


256 x 256 


B1S9 


1 


0.9 


- 


- 


- 


OFF 


256 x 256 


B4S0T3M9C 


4 





3 


10 -9 


2.88±2.02 x 10 -9 


ON 


256 x 256 


B4S5T3M9C 


4 


0.5 


3 


10" 9 


5.33±4.87 x 10~ 9 


ON 


256 x 256 


B4S7T3M9C 


4 


0.7 


3 


10" 9 


5.38±4.06 x 10 -9 


ON 


256 x 256 


B4S98T3M9C 


4 


0.98 


3 


10" 9 


3.98±3.36 x 10~ 9 


ON 


256 x 256 


B4S9rT3M9C 


4 


-0.9 


3 


10 -9 


0.64±0.47 x 10 -9 


ON 


256 x 256 


B4S9T1M9C 


4 


0.9 


1 


10 -9 


3.85±3.50 x lO" 9 


ON 


256 x 256 


B4S9T10M9C 


4 


0.9 


10 


10" 9 


3.62±3.90 x 10 -9 


ON 


256 x 256 


B4S91 


1 


0.9 


- 


- 


- 


OFF 


192 x 128 


B4S9h 


4 


0.9 








OFF 


384 x 384 


B4S0T3M8C 


4 





3 


10 -8 


3.73 ±2.18 x 10~ 8 


ON 


256 x 256 


B4S9T3M8C 


4 


0.9 


3 


lO" 8 


4.99 ±3.80 x 10~ 8 


ON 


256 x 256 


B4S9T3M7C 


4 


0.9 


3 


6.3 x 10" 8 


1.69 ±1.39 x 10~ 7 


ON 


256 x 256 


B4S0T3M7C 


4 





3 


6.3 x 10" 8 


2.16 ±1.46 x 10~ 7 


ON 


256 x 256 


B4S5T3M7C 


4 


0.5 


3 


6.3 x 10" 8 


2.03±1.85 x 10~ 7 


ON 


256 x 256 


B4S75T3M7C 


4 


0.75 


3 


6.3 x 10" 8 


1.65±1.15 x 10~ 7 


ON 


256 x 256 


B4S98T3M7C 


4 


0.98 


3 


6.3 x 10" 8 


1.37±1.08 x 10~ 7 


ON 


256 x 256 


B4S9rT3M7C 


4 


-0.9 


3 


6.3 x 10" 8 


0.39±0.41 x 10~ 7 


ON 


256 x 256 


B4S9T1M7C 


4 


0.9 


1 


6.3 x 10~ 8 


1.53±1.02 x 10~ 7 


ON 


256 x 256 


B4S9T10M7C 


1 


0.9 


10 


6.3 x lO" 8 


1.43±0.88 x 10~ 7 


ON 


256 x 256 



a Disk scale height of H/r w 0.07 instead of 0.12 as for B4S9T3M9C 

^ Anag.o = 50 instead of 10 

c r max = 1.1 X 10 4 r<3 instead of 120rQ 
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Figure 5. Broadband spectra comparing, at three different mass 
accretion rates: 10 — 9 MQ/yr (black), 10 -8 M0/yr (blue), and 
10~ 7 M0/yr (orange), the effect of turning ON (solid) and OFF 
(dash-dot-dot) the cooling function in simulations with the same 
set of initial parameters (4-loop model, a* = 0.9, Ti/T e = 3 and 
i = 85°). Emission from simulations B4S9T3M8 and B4S9T3M7 
have been respectively multiplied by 10 and 1000 to improve read- 
ability. 



nitudes) appear at a mass accretion rate of 10 8 Mq/jt and 
M = 1(T 7 M Q /yr. 

Next, we have naturally focused our attention on how 



Figure 6. Broadband spectra comparing the effect on emission 
of varying the mass accretion rate, when the cooling function is 
turned on . All simulations have been set with an initial 4-loop 
magnetic field, a black hole spin a* = 0.9, a temperature ratio 
Ti/T e =3 and an inclination angle i = 85°. 



varying the mass accretion rate changes the resulting spec- 
tra. As shown in Figure [6] which compares SEDs from simu- 
lations B4S9T3M9C, B4S9T3M8C and B4S9T3M7C, there 
is a significant positive correlation between the mass accre- 
tion rate and the emission at all wavelengths. 
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Figure 7. Broadband spectra of simulations with the same set 
of initial parameters (4- loop model, Tj/T e = 3, mass accretion 
rate M = 10 -9 M Q /yr and i = 85°) but for different black hole 
spin values: a* = {0.7,0.9,0.98} (top panel), o» = {-0.9,0,0.5} 
(bottom panel). 



Figure 8. Broadband spectra comparing the effect of varying the 
ion-to-electron temperature ratio for simulations with an initial 
4-loop magnetic field, a black hole spin a* = 0.9, a mass accretion 
rate M = 10 — 9 A4"Qyr~ 1 and an inclination angle % = 85°. 
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We also test the role of black hole spin a,. As can be 
seen from Figure [7] the luminosity rises from the case of a 
retrograde spinning black hole to a spinning one of a* = 
0.5. Then, it appears that the emission reaches a plateau 
between a* = 0.5 and a, = 0.7, before rising again. More- 
over, simulations from models with spin a» = -0.9, 0, 0.5 and 
0.98 seem to be more variable compared to models of spin 
a* = 0.7 and 0.9, which lead to a higher variability in the re- 
sulting radiation. Finally comparing emission from positive 
and negative spins reveals almost four orders of magnitude 
difference between fluxes at same absolute spin value. 

Regarding the ion-to-electron temperature ratio, it is 
straightforward to assess how this parameter affects the ra- 
diation in non-cooling simulations, but it is not for cooling 
ones. On one hand, in non-cooling simulations, the temper- 
ature of the ions stays the same. So increasing the temper- 
ature ratio decreases the temperature of the electrons and 
therefore the emission. On the other hand, in cooling simula- 
tions, effects from radiative losses and from efficiency of cool- 
ing processes between ions and electrons conflict and their 
results on spectra are not straightforward. Radiative losses 
will lower the temperature of the electrons while increas- 
ing Ti/T e results in having less efficient cooling processes, 
therefore less radiative emission from the electrons and thus 



Figure 9. Broadband spectra comparing the effect of varying the 
initial magnetic field configuration on the emission of simulations 
with a black hole spin a» = 0.9, a temperature ratio Ti/T e = 3, a 
mass accretion rate M = 10 — 9 M0yr — 1 and an inclination angle 
i = 85°. 

a higher electron's temperature. Nonetheless Figure [8] shows 
that in our specific case, increasing Ti/T e decreases the total 
emission. However, at higher accretion rates this increase is 
less than the F„ ~ M 2 scaling without cooling. 

Our work is one of the first to address the question of 
how magnetic field configuration model in the initial accre- 
tion disc affects the resulting emission, in the case of Sgr A*. 
Figure [5] compares spectra from two magnetic field configu- 
rations: the 1-loop and the 4-loop models. The figure shows 
that, while in the submm and the near-infrared bands the 
emission is fairly independent of the model - both spectra 
are within each others variability range - in the X-ray, the 
emission is very sensitive to the initial configuration. 

The inclination angle is the last parameter we tested. 
This parameter is only used in the post-processing ray- 
tracing code grtrans. It defines the viewing angle of a 
distant observer on the system. Doppler shifts and optical 
depth are the two signifi cant changes induc ed by varying 
the inclination angle (e.g. jDexter et al]|2009l ). For most ob- 
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Figure 10. Broadband spectra comparing the effect of varying 
the viewing inclination angle on the emission of the reference 
simulation B4S9T3M9C. 





B4S9T3M9C 

B4S7T3M9C 








1/ 

I / 




1 
1 


10 10 10" 


10 12 10 13 10 14 
v[Hz] 


10 15 



Figure 11. Broadband spectra of our preferred parameters fits. 
The models suggest that the mass accretion rate at which Sgr A* 
accretes is around 10 _9 MQ/yr and the spin of the central black 
hole is likely to be between a* = 0.7 and a* = 0.9. 



4.06 x 10~ 9 A/ Q /yr in simu lation B4S7T3M9C. These values 
confirm the conclusion of |D12| that it is not necessary to 
consider radiative processes when simulating accretion onto 
Sgr A* because of its exceptionally low accretion rate. 

Our favored target mass accretion rate is consistent 
with the lower limit imposed by observations of linear po- 
larisation in the submm bump. It is interesting to note that, 
at this mass accretion rate, two models provide a good de- 
scription of the data, with their only difference being the 
spin parameter value. Figure [TT] shows that our models with 
a, = 0.7 and a* = 0.9 match Sgr A* data at 230 GHz and 
~ 5.6 x 10 34 erg/s. These two models distinguish themselves 
only in the description of the near-infrared observations. 
Our results are consistent with the upper limit of a* = 0.86 
at 2cr significance given bv lBroderick et all (1201 if) obtained 
with millimetre - VLBI observations |Doeleman et al.l 120081 : 
iFish et alJlioTll ). 

As reported in lDexter et alJ l)2010h . the favored ion-to- 
electron temperature ratio depends strongly on initial con- 
ditions since the temperature of the ions scales with the 
disc thickness. The real constraint is therefore on the tem- 
perature of the electrons. Currently, our ability to constrain 
this parameter is only as good as the code itself. Although 
our preferred parameter-space fits are found with a ion-to- 
electron temperature ratio of 3 which suggests that the pro- 
cesses coupling the ions to the electrons in the accretion 
disc are mildly inefficient, we would advise caution regard- 
ing this conclusion. Similarly, no strong constraints can be 
drawn from our study of the inclination angles. Our pre- 
ferred parameters fit is obtained for an inclination angle of 
85 ± 40 deg, which is not unexpected given our position in 
the plane of the Galaxy. 

A further interesting point is that the X-ray emission 
is very sensitive to the initial magnetic field configuration 
in our simulations, while the submm emission is fairly in- 
dependent. The X-ray upper limit of Sgr A*'s data may be 
a promising way to constrain the effect of magnetic fields. 
We only tested a limited set of initial conditions, and in this 
context, we obtained our best fit when seeding the initial 
torus with a four sets of poloidal magnetic field loops. 



server inclinations, Doppler beaming is the predominant ef- 
fect from ray-tracing from Keplerian discs while at lower in- 
clinations optical depth is the main effect. Figure [lOl shows 
that, at higher inclinations, Doppler beaming leads to larger 
fluxes, and moves the peak of the spectrum to higher fre- 
quency. At lower inclinations, the figure also shows that op- 
tical depth causes larger variability of the overall fluxes. 

5 DISCUSSION 

5.1 Preferred parameter space 

The present study was designed to determine which set of 
parameters will be the closest to reproducing the quiescent 
state of Sgr A*. Our results show that the most compatible 
spectra with the observational data are those of simulations 
B4S9T3M9C and B4S7T3M9C. We find that we can fit the 
Sgr A* data at 230 GHz, which has a value of 3 Jy, or 
vL v ~ 5.6 x 10 34 erg/s, with a mass accretion rate of 2.60 ± 
1.54 x 10~ 9 Af Q /yr in simulation B4S9T3M9C and 5.38 ± 



5.2 Cases of retrograde spin 

Figure [12] shows the spectra of retrograde spin models at 
mass accretion rates of 1O _9 M0 yr _1 and lO _7 M0yr _1 . It 
is interesting to note that the spectral shapes of these mod- 
els are significantly different from the positive spin models. 
The most notable difference is that the synchrotron emis- 
sion peaks at higher frequency. This peak originates from 
an emitting region inside 3r g , filled with material at a very 
high temperature (~ 10 13 K). While this very high temper- 
ature and the displacement of the synchrotron peak might 
be seen as criteria to exclude a* < models for Sgr A*, 
based on the submm VLBI data, it is likely premature to 
draw any definite conclusions. Further analysis is required 
to understand the origin of the very hot material and to es- 
tablish whether a retrograde spin model for Sgr A* with a 
high mass accretion rate is viable. 
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Figure 12. Broadband spectra of retrograde spin simulations 
with the same set of initial parameters (4-loop model, Ti/T e = 3 
and i = 85°) but at different mass accretion rates. 
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Figure 13. Light curves at three different bands: 2.1 X 10 17 Hz 
(top), 4.5 X 10 14 Hz (middle) and 9.4 X 10 11 Hz (bottom) of the 
simulation B4S9T3M9C. 



5.3 Flaring events 

Although the main goal of our work was to reproduce the 
quiescent state of Sgr A*, we want to say a few words about 
the origin of flaring events leading to fast variability within 
our spectra. As shown in Figure 1131 our preferred simula- 
tion (B4S9T3M9C) experiences a brief X-ray flare, with the 
emission in this band increasing by three orders of magni- 
tude over ~ 20 minutes. It is apparent from Figure [T?] which 
presents the radial and the 6 angle profile of the peak of the 
X-ray emission, that this emission originates from a narrow 
region located between 2.3 rad < 9 < 2.8 rad, composed of 
two blobs: the first situated around r — 5 r g and the second 
spreading between 10 r g and 14 r g . 

To investigate the origin of the blobs, we examined 
temperature maps at each time step of the simulation to- 
gether with the evolution of its magnetic field lines. Fig- 
ure [15] presents four snapshots of the formation of the blobs. 
The initial state of the flare region is a thin filament of co- 
herent field starting at the event horizon and extending out 
to r ~ 15 r g , which sug gests that the simulation develops a 
channel-mode solution (jrlawlev fc Balbuslll992D in this re- 
gion. 

These episo dic flares in our 2.5D sim ulations are very 
similar to what iDodds-Eden et all (|2010t) report in their 
study of large, sporadic magnetic reconnection events oc- 
curring near the last stable circular orbits in their 2.5D 
GRMHD simulations. They suggest that because these 
events have timescales and energetics consistent with Sgr 
A*'s flares, they may represent actual physical mecha- 
nisms. However, Sgr A*'s X-ray flares always have a si- 
multaneous infrared counter-part, while only the largest in- 
frared flares show X-ray flares in general l|Eckart et al.|[200fj : 
IDodds-Eden et al"1l201ll ). Our simulated light curves do not 
show an infrared event corresponding to the X-ray flare. 
Moreover, no similar behaviour has ever been reported in 
3D GRMHD simulations thus it is highly likely that these 
flaring events are numerical artefacts rising from the two di- 
mensional nature of the simulations when magnetic recon- 
nections occur near the event horizon of the black hole. By 
enforcing axisymmetry, 2.5D simulations allow larger coher- 
ent magnetic field structures to form, enhancing variability 




r(M) 




1.0 1.5 2.0 
e (radians) 



Figure 14. Radial and profile of the logarithm of the X-ray 
emissivity, area and redshift weighted. The emissivity is calcu- 
lated at the peak of the X-ray flares shown in the top panel of 
Figure [13] 



when these structures finally reconnect. We choose to min- 
imize the effect of these rare events on the final spectra by 
using median rather than time-averaged spectra as discussed 
above. 



5.4 Comparison with previous works 

In the past few years, several groups have focused their anal- 
ysis on comparing simulated observation of Sgr A* to data. 
While the general setup of our simulations are similar to 
that of other groups, the details of the code itself and the 
treatment of radiation are not. Somewhat surprisingly, and 
encouragingly, this study produces results which corrobo- 
rate the findings of a great deal of the previous work in 
this field, suggesting that all groups are converging on a 
consistent picture for conditions around the supermassive 
black hole. In particular, our work is in agreement with the 
findings of iMoscibrodzka et al.l (|2009l ) which showed that a 
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Figure 15. Snapshots of a flaring event occurring in simulation 
B4S9T3M9C at 11.73h (top-left), 12.381h (top-right), 12.754h 
(bottom-left) and 12.847h (bottom-right). Each snapshot shows 
a map of the temperature (color) and magnetic field lines (black 
lines) . 



radiative model of Sgr A* with a mass accretion rate of 
1.86 x 1O _9 M0 yr _1 , an ion-to-electron temperature ratio of 
3, a spin of a, = 0.94 and an inclination angle of 85 ± 40 deg 
fits Sgr A* observations. 

These findings are also consistent with those of 
iDexter et all (|200sl . l2010l ) and further support the idea that 
the submillimeter bump in Sgr A*'s data originates from 
within the innermost region of an accretion disc, accreting 
on Sgr A* at a mass accretion rate of ~ 2 x 10 _9 Mq yr -1 . 
However, it is worth noting that all simulations so far are 
finding the submm bump to be dominated by the emission 
from the inner disc because the jets are not yet correctly 
physically described. Most importantly, idealized MHD pre- 
vents realistic mass loading in the jet funnels, but the res- 
olution of the grid is often poor along the poles, and the 
numerical floors often dominate in these regions of the sim- 
ulations as well. 

Another important result from this work is a quantita- 
tive measure of the increasing importance of a self-consistent 
treatment for radiative cooling losses in GRMHD simula- 
tions, with increasing mass accret ion rates. This result sup- 
ports our companion paper |P12r s conclusion that above a 
mass accretion rate of ~ 10~ 7 A/Edd, a self-consistent treat- 
ment of the radiative losses in GRMHD simulations not only 
affects the dynamics of the simulations, it also affects the 
radiative emission. This conclusion will affect any previous 
works done on Sgr A* which neglected the radiative losses 



and used a mass accretion rate hi gher than this limit. For 
example, IShcherbakov et all l|2010l ) fits Sgr A* data in the 
submillimeter bump and uses polarised radiation to find a 
mass accretion rate of (1.4 — 7.0) x lO _8 M0yr _1 . Our result 
here implies that their final spectrum would be affected by 
cooling losses. 



5.5 Limitations 

The current study has significant limitations. All of the 
GRMHD simulations presented here are axisymmetric 
(2.5D). Axisymmetric simulations cannot sustain turbulence 
and so never reach a quasi-steady state. Axisymmetry also 
tends to exaggerate variability relative to the 3D case, and is 
likely responsible for the rare, large amplitude flaring events 
seen in many of our simulations. 

Another limitation of our study, shared in general by the 
current class of ideal MHD simulations, is that the jets can- 
not be mass-loaded. Observations of flat/inverted spectra 
from compact jets in LLAGN, indicate optical depth effects 
which the current simulations cannot approach. Most likely 
once prescriptions for mass loading and particle acceleration 
in the jets are included, emission from the base of the jets 
will increase in the submm for Sgr A*, and have some effect 
on our favored parameter space. 

We have also simulated a limited set of initial condi- 
tions. We have found that the initial magnetic field configu- 
ration can have an important effect on the resulting spectra, 
especially at high energies, but a wider range of configura- 
tions should be t r ied to fully explore this issue. For instance, 
iMcKinnev et all (|2012l ) argue that the initial condition used 
in these simulations artificially restricts the available mag- 
netic flux, and show that large amounts of coherent flux can 
significantly alter the dynamics of the accretion flow. 

Finally, this and almost all previous studies attempting 
to constrain the parameters of Sgr A* have assumed that the 
accretion flow angular momentum axis is aligned with the 
black hole spin axis. However, th i s is un likely to be the case 
in reality, and IDexter fc Fragile! ([2012) show that spectral 
fits can change dramatically even for tilts as small as 15°. 



6 SUMMARY 

This paper presents for the first time self-consistent spectra 
from radiatively cooled GRMHD simulations of the accre- 
tion flow around a black hole, in particular, Sgr A*. Our 
study concludes that the central black hole is most likely 
rapidly spinning (0.7 < a* < 0.9), and that Sgr A* is accret- 
ing at a mass accretion rate of ~ 2 x lO _9 M0yr _1 . While 
no significant conclusions can be drawn from the correla- 
tion between the resulting emission and the initial magnetic 
field configuration model, we obtain our best description for 
the submillimeter data by seeding the initial torus with a 
4-loop poloidal magnetic field, suggesting that a more com- 
plex morphology could be favoured. Finally no constraints 
on the inclination angle can be derived from our work, but 
it is consistent with the general sense that Sgr A* should be 
more edge on than face on. 

Our work confirms the limit on the mass ac cretio n rate 
(~ 10 -7 Msdd) reported in our companion paper If5l2l . where 
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self-consistent treatment of cooling losses in GRMHD sim- 
ulations becomes important. Above this limit, spectra gen- 
erated from GRMHD simulations where radiative losses are 
not taken into account can be potentially orders of magni- 
tude too high. However, for other sources the exact limit 
may vary slightly with the mass and spin as well as initial 
conditions of the simulation. Nonetheless, this result is very 
important to keep in mind for future studies of more typical 
nearby LLAGN such as M81, M87, etc.. 

We showed that high energy emission from GRMHD 
simulations is sensitive to the magnetic field configuration 
in the initial accretion disc. Further research regarding the 
role of the magnetic field configuration in the dynamics and 
radiation of GRMHD simulations may ultimately help dis- 
tinguish between models for the origin of the magnetic fields 
close to the black hole. 

Recently there have been claims o f a 3 pc-scale, jet- 
drive n outflow from Sgr A* in the radio (|Yusef-Zadeh et aU 
|2012| ) as well a large-scale jet fe ature in the Fermi GeV 7- 
ray maps of the Galactic center l|Su fc Finkbeine"rll2012t ). If 
one or both can be confirmed, these features will provide 
valuable constraints for the next technological development 
of GRMHD simulations, which is the inclusion of more real- 
istic mass-loading and particle acceler ation in the jets. Sim - 
ilarly, the discovery of the G2 cloud (|Gillessen et alj|2012h 
on a collision course with Sgr A* for 2013 may provide new 
tests of Sgr A*'s emission at higher accretion rates, for com- 
parison with our results here. 
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